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We review and refine the concept of a mean-fieid theory for the study of sandpile models, which 
are of central importance in the study of self-organized criticality. By considering the simple one- 
dimensional random walker with an absorbing and reflecting boundary we are able to construct a 
complete mean-field picture which we can solve in detail for different types of driving. Using this 
theory, we are able to clarify the effect of finite driving rate on sandpile models, as well as the 
observed sensitivity of certain universal quantities on the driving. 
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Sandpile models form the core of the study 
of self-organized criticality (S0C)[l-4]. These 
models describe the movement of particles in 
bursts called avalanches, with SOC character- 
ized by the fact that the avalanche size prob- 
ability exhibits simple finite-size scaling, which 
can be described by the critical exponents r 
and A (defined below). These exponents have 
been determined exactly for very few models 
[2], with the exceptions generally being mean- 
field theories, where it is well established that 
T — 2>/2 and A = 2 [5, 6]. However, one ques- 
tion which has not been addressed in the lit- 
erature is whether mean-field theory can tell 
us anything more qualitative about the nature 
of sandpile models. In the following we ad- 
dress this question by formulating and solving 
a mean-field theory based on escape times of a 
random walker. We find, of course, that there 
is scaling and universality in this model, and 
that it arises as a direct consequence of symme- 
try and bulk conservation. We also see clearly 
the roles played by the exponents r and A and 
clarify the sensitivity of these exponents on the 
driving, as observed numerically in Refs. [4, 7- 
11]. Finally, we identify the effect of non-zero 
driving rates, which has been of recent interest 
because of connections to absorbing state phase 
transitions [11, 12], and use this as a basis for 
writing down the scaling theory of SOC. 

The difficulty with solving sandpile models 
analytically is that there is a lattice which can 
contain particles which are stable as well as 
unstable. During the course of an avalanche, 
only unstable particles move but stable parti- 
cles may be made unstable by interacting with 
those which are already unstable. Thus, the 
stable particles form a 'medium' which interacts 



with the avalanches formed by unstable parti- 
cles. Keeping track of these interactions soon 
becomes analytically intractable and thus there 
is a real deficiency in exact solutions for sand- 
pile models. A mean-field theory replaces the 
medium seen by each degree of freedom with 
a self-consistent average, such that fluctuations 
are ignored. In this case, we look at a moving 
particle as it enters a site. There will be some 
probability that it is made stable, as well as a 
probability that it excites a stable particle al- 
ready there to become unstable. If the sandpile 
is in a stationary state and particles are con- 
served in the bulk, then the average number 
of unstable particles that leave a (bulk) site is 
equal to the average that entered. Hence, in the 
mean-field theory, we neglect the interactions 
of the particles entirely and unstable particles 
simply perform a random walk until they leave 
the system through a boundary site. 

To define our mean-field theory, we consider 
a random walker on a line, x G [0, L] with an 
absorbing boundary a,t x — Q and a reflect- 
ing boundary a,t x = L. We introduce the 
walker at some point < xq < L and define 
the 'avalanche' size as the time, T, it takes for 
the walker to leave the system through the ab- 
sorbing boundary at a; = 0. The probability 
density for the walker at time t, which starts at 
Xq at t = is (f) = 4>{x, t; xq) and the dynamics 
are described by the diffusion equation 



d(j) £)^^'t' 
dt dx"^ 



dx ' 



(1) 



where D is the diffusion constant and v is the 
drift velocity. The probability density for T is 
simply given by the current out of the system 
at time t — T: 



P{T-L) = D 



dcj){x,T]XQ) 
dx 



(2) 
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As well as being a mean-field theory for sand- 
pile models like those in Refs. [1-4], this model 
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can be justified as a mean-field theory of the 
Eulerian walker [13] and is in the same uni- 
versality class as the critical branching process 
where each tree is restricted to a maximum of 
L generations [14. 15]. 

The SOC hypothesis predicts that under cer- 
tain circumstances the event-size probability 
exhibits simple finite-size scaling, that is, for 
T » 1 and L > 1, 



P{T- L) = aT-'"S{hT/L^), 



(3) 



where a and b arc non- universal constants, r 
and A are universal exponents and is a scal- 
ing function. We can also study the moments 
of this probability density 

poo 

(T") = / T"P{T] L)dT (4) 
Jo 

which will scale with system size 



(T") « r„L^" for L > 1, 



(5) 



where r„ are non-universal amplitudes and 

7„ — A(l + n — t). For the mean-field the- 
ory we will identify the reason for this scaling, 
which may give us clues to the cause of scal- 
ing in SOC systems. We shall also yield some 
explanation for why in SOC systems wc find 
that T depends on the location of the driving, 
Xq, whereas A remains constant: Stationarity 
gives us the scaling relations (for instance, see 
Ref. [7]) 

A(2-r) = ^ -o-'^o^^t- (6) 
^ ' \2 xo/L = const. ^ ' 

and it is found numerically that A is the same 
for both cases whereas r changes [4, 7-11]. This 
is important to understand as we might ex- 
pect from the renormalization group that uni- 
versal quantities such as critical exponents are 
only affected by symmetry-breaking or non- 
conservative perturbations. 

For now we set the drift velocity w = 0, an- 
ticipating the fact that scaling is observed only 
in this case. We shall handle the finite drift 
case later in this letter. If the diffusers start 
close to the absorbing boundary then they take 
a time « L^/D to reach the reflecting bound- 
ary on the right hand side and for small times 
the system does not 'feel' that the system is fi- 
nite. Hence, when t <C Lt^/D, (l){x,t;xo) takes 
the same form as the propagator for diffusion 
on the semi-infinite line. 
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(x-xq) 

" — 4bt 



e 



(7) 



For t ^ L'^/D, only the lowest order mode of 
^{x,t;xo) survives and 



2L3 



xoxe 



-Dt 



(8) 



Taking the derivative in Eq. (2), we find 

,T-i xl/D « T < L'^/D 

T > L'^/D 



P{T-L) 



which can be put in the scaling form 
P{T; L) « xoT-i<^{DT/L'^) 



(9) 



(10) 



with ^^(x) = const, for small arguments and 
2.3/2g-2: large arguments. For constant xq, 
this is simple finite-size scaling and by compar- 
ing to Eq. (3) we identify A = 2 and t = 3/2. 
In the case of bulk driving, xq = yoL, where yo 
is a constant, we have to recast the equation 



P{T;L) w yoLT-^''^<^{DT/L^) 
= yoT-^nDT/L'') 



(11) 



where ^{x) = x~^^'^^(x). Hence, for oc i, 
we still have A = 2 but r = 1. 

We also find the scaling behavior of the mo- 
ments. This can be done by writing down the 
generating function, Z = Z{X,xo; L), where 



(T") = (-1)" 



A=0 



This is easily determined to be 

j^2n sinh\/Acosh ^(1 - ^)VX 



Z = 2- 



sinh2\/A 



(12) 



(13) 



We note that (T") has dimensions of (time)" 
and dimensional analysis tells us that 

r2n 

(T") = —gn{xo/L) (14) 

where <?„ is a scaling function for the nth mo- 
ment of the dimensionless ratio xo/L < 1. If 
we write down <;„ as a power series in xo/L, we 
have for n > 0, 



(15) 



where the coefficients C^.n are numbers inde- 
pendent of L, D and xq- It is clear that the 
leading order behavior in L comes from the low- 
est order term, so that for n > 0, 



r2n-l 



(16) 
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We can find the Ci,„ directly from Z: 



(2n)! 



(17) 



where -B„ are Bernoulli numbers. 

Prom Eq. (15) we can write down the scaling 
of the moments in the usual scaling form, 



(18) 



with 7„ = A(l + n — r) where A = 2 is the 
gap exponent and r = 3/2 is the escape time 
exponent. The amplitudes r„ are non-universal 
and are equal to 



(-l)"-^n! 



22"(2^" - \)D2 
(2n)! 



(19) 



Note that for bulk driving x^jL = const., 
Eq. (16) implies that t changes its value to 1, 
but A remains unchanged, in agreement with 
Eq. (11). 

We also consider the moment ratios 



an 



{T'^){Ty 

(T2)"- 



rj^n— i 
1 

-pn— 1 



(20) 



which are defined for n > 0. For an event-size 
probability given by Eq. (3), the Qn approach 
universal values for L ::J> 1. However, we see 
in Eq. (15) that for xq/L — const., each term 
in the expansion of (T") will be of the same 
order in L, and so the amplitude r„ will be the 
sum of all the terms and not just the first. We 
therefore expect universal amplitude ratios to 
differ for the two types of driving. 

We can therefore see that the reason for the 
dependence on r and not A on the driving is 
because the probability density, P{T;L), and 
therefore its moments are proportional to xq 
to lowest order. In terms of the exponents 7„, 
this means that the gap between the 7n, which 
is proportional to A, is unchanged. We argue 
that the linearity in xo can be shown to be a 
consequence of conservation of particles, as well 
as a symmetry between x and xq in diffusion: 
(f){x, t; xo) not only obeys the diffusion equation 
written in Eq. (1), but also 



d(j) B'^cj) d(j) 



(21) 



One interpretation of this is as a time reversal 
symmetry: the probability of a particle diffus- 
ing from a; to y in time t is identical to that 
for y to X, providing we reverse the direction 
of the drift. Even the presence of absorbing 
boundary conditions does not break this sym- 
metry. The effect of this is clear. Symmetry 



and conservation constrain the dependence of 
the observables on xo to be the same as for x, 
which is linear to lowest order. 

There has been much recent interest in study- 
ing the effects of finite drive rate on the behav- 
ior of sandpile models. It is widely held that 
scaling only occurs when the driving rate is set 
to zero, with a finite drive rate introducing a 
cutoff length scale for avalanche sizes. In our 
mean-field theory, we wait until each walker 
leaves the system before adding the next, which 
corresponds to an infinitely slow driving rate. 
Consider now the situation where we add walk- 
ers at some finite rate q and measure the event 
sizes as the time between instances of having 
zero walkers in the system. This situation is 
modeled well by treating q like a drift, q ^ —v. 
Taking the solution in the limit L — > oo, we find 
for T > xl/D, 



P{T;q). 



Xo 



3 

■26 



(22) 



Hence, there is an exponential cutoff with the 
characteristic time 



rr, 4£» 

T = 

r 



(23) 



and the model exhibits scaling only when Tq 
oo, that is, q 0. This is in agreement with 
the idea that SOC is only achieved in the limit 
of zero drive rate [12]. 

We are now in a position to formulate a gen- 
eral scaling theory for the moments of the event 
size probability in sandpile models. We have a 
number of length scales, L, xq and D/q, which 
we have already introduced, plus an additional 
length scale corresponding to details of fiuctu- 
ations on the lattice, a. Writing (T") in terms 
of dimensionless quantities, we have 



L' L' D 



(24) 



where is some scaling function. The scaling 
function may scale with a, but we do not expect 
it to produce a cutoff in the scaling function. 

Hence, we write 

where r]{n) is an exponent which may depend 

on n and S„ is the scaling function $„ with the 
a dependence factored out. For sandpiles with 
conservation we anticipate similar constraints 
on Xo as for the mean-field model and hence we 
make the conjecture 

(T") = xoL^^-^'-^^-'^^n (§) (26) 
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where, again, is obtained from S„. It is now 
clear that when q ^ 0, we lose another length 
scale and (T") exhibit pure scaling. 

We have investigated the mean-field theory of 
sandpiles in some detail in order to gain insight 
into the nature of the "criticality" of SOC sand- 
piles. We have found that a simple mean- field 
theory exhibits scaling in the presence of bulk 
conservation and slow driving. We have shown 
that critical exponents and amplitude ratios, 
are dependent on the type of driving used, a 



feature which has been observed numerically in 
non-trivial sandpile models. It should be noted 
that this mean-field description depends only 
on conservation and symmetry in the underly- 
ing sandpile model, and does not describe any 
mechanism of self-organization or fluctuations 
in the medium, which are the unique proper- 
ties of sandpiles. However, it is remarkable that 
such a simple mean-field theory has been able to 
capture so successfully the remaining features 
of these non- trivial models. 
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